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Abstract 

Topological invariants for the AD gapped system are discussed with 



c 



O ; application to the quantum vacua of relativistic quantum fields. Ex- 



pression A/3 for the AD systems with mass gap defined in [13] is con- 



D . sidered. It is demonstrated that A/3 remains the topological invariant 



when the interacting theory in deep ultraviolet is effectively massless. 
We also consider the 5D systems and demonstrate how AD invariants 
emerge as a result of the dimensional reduction. In particular, the 
new AD invariant A/5 is suggested. The index theorem is proved that 
defines the number of massless fermions in the intermediate vac- 
uum, which exists at the transition line between the massive vacua 
with different values of A/5. Namely, 2np is equal to the jump AA/5 
across the transition. The jump AA/3 at the transition determines the 
number of only those massless fermions, which live near the hypersur- 
face oj = 0. The considered invariants are calculated for the lattice 
model with Wilson fermions. 



1 Introduction 



The momentum space topology is becoming the main tool for the investiga- 
tion of the robust properties of ground states (vacua) of condensed matter 
systems (see review papers [U EJ [3]). It revealed variety of vacuum struc- 
tures with nontrivial topological invariants, which is now commonly known 
as topological materials. These invariants protect gapless fermions in bulk, or 
on the surface of the fully gapped topological materials, or inside the core of 
topological defects in these systems, such as strings, vortices, monopoles, do- 
main walls, solitons, etc. The gaplessness (masslessness) of fermionic species 
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is not sensitive to the details of the microscopic physics: irrespective of the 
deformation of the parameters of the microscopic theory, the energy spec- 
trum of fermions remains strictly gapless. Though the physics of topological 
materials started 30 years ago with the pioneering work by Nielsen and Ni- 
nomiya [I], who demonstrated that the masslessness of elementary particles 
can be protected by the momentum space topology, the investigation of the 
topologically nontrivial vacua in relativistic quantum field theories is still 
in its infancy, with a rather few papers exploiting topological invariants in 
momentum space (see e.g. [HI El [Tl El [HI [101 [111 H2] ) - 

Here we are interested in the topologically nontrivial relativistic vacua 
with massive fermions. The AD vacua with mass gap may exist in vari- 
ous phases, which are described by different values of topological invariants. 
Recently such an invariant was considered in [13]. This invariant A3 is deter- 
mined on the to = hypersurface of AD momentum space, which represents 
the 3D Brillouin zone. Such invariant is not sufficient for the full classification 
of the topological states of the lattice models, whose momentum-frequency 
(a;, p) space represents the AD Brillouin zone. Here we introduce another 
topological invariant A/5, which is appropriate for the lattice models (see 
Eq.( )20|) ). The invariants A/3 and A/5 are then applied to particular lattice 
model with Wilson fermions (see e.g. Ref. [H]). Table |2] demonstrates the 
values of these invariants in different ranges of mass parameter m, in which 
the vacuum represents different phases of topological insulator. 

There are crtitical values of mass parameter in the Wilson model m = 0, 
—2, —4, —6 and —8, at which the topological quantum phase transitions oc- 
cur between the insulators with different A/3 or/and A/5 (more on topological 
phase transitions, at which the topological charge of the vacuum changes 
while the symmetry does not, see [35]). At these values of m the vacuum 
states are gapless and they represent analogues of topological semimetals: 
similar to the four- dimensional graphene discussed in Ref. [11] they contain 
massless Dirac fermions. There is an analogue of index theorem, which re- 
lates the number of these massless fermions with the jump of the topological 
invariant across the transition (AA5 or AA/3). The total algebraic number 
of gapless fermions at the transition point is np = AA/5/2, and according 
to Table [2] the states with critical m — 0, —2, —4, —6 and —8 have corre- 
spondingly 1, 4, 6, 4, 1 massless species. The jump AA/3 also determines 
the number of massless fermions in the intermediate states, but only those of 
them which live near the hypersurface 00 = 0, i.e. np(u = 0) = AA/3/2 [13J. 
The intermediate states have correspondingly 1, 3, 3, 1, of such massless 
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species. 

Since the invariants are expressed in terms of the Green's function, they 
are applicable to the interacting systems as well [9]. We prove that the 
relation between the number of gapless fermions and the jump in the topo- 
logical charge remains valid within the interacting theory for the wide class of 
4D models. (Recent discussion on topological invariants for interacting con- 
densed matter systems in terms of Green's function and some peculiarities 
in the bulk-boundary correspondence there see in Refs. [T6l ITT] .) 

The paper is organized as follows. In Section [2] we consider topological 
invariant A3. In Section [3] we consider the 5D massless systems. In Section 
H] we derive topological invariants in 4D via the dimensional reduction of 
the 5D constructions. In Section we suggest the method of calculation of 
the invariants A/5 and A/3 and formulate the index theorem. In Section [6] we 
apply the suggested constructions to the lattice model with Wilson fermions. 
In Section [7] we end with the conclusions. 

2 Topological invariant A/3 
2.1 Green functions 

Let us consider for the defmiteness K flavors of Dirac fermions x A ,A = 
1, K coupled to some gauge field A, that, in turn, may be coupled to the 
other unknown fields. Most of the conclusions, however, may be extended to 
a more general case. The fermion Green function in Euclidean space has the 
form: 

= \ J D x DAexp(-S G [A] 

- J d 4 x X m + iAiW + m\x)xA(0)x B (x) 

= 1 J DxDAexp(-S G [A}) Det + iAi)Y + H 
{m+iAiW+im]- 1 ^ (1) 

Here Sq is some unknown effective action. We use Euclidean 7 - matrices 
such that {7*,7' J } = 25 1 ^ . It will also be implied that a certain Lorentz 
invariant gauge is fixed. When the interactions are turned off, we have 

G(u, p)5 = / d 4 *£f (0, x)e^ 4+ ^ = ^ • » ^ = 00 (2) 
J p,x — im 
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When the interactions are turned on, we suppose that Lorentz invariance 
and parity are not broken. We shall also use the function Q = —i^Q that 
is Hermitian as follows from ([1]). It can be diagonalized with respect to the 
flavor index via Unitary transformations of the fermion field. Later we shall 
imply that the basis is chosen, in which Q (and Q) is diagonal in flavor index. 
Symmetry fixes the form of the Green function 

G(u,p)% = Z c A (p 2 )[ . 1 . lg (3) 

Here Z(p 2 ) is the wave function renormalization function, while m(p 2 ) is 
the effective mass term. Both Z and m are diagonal matrices in indices A, B. 
Below we omit fermion indices A, B in most of the expressions. 

In addition we shall also consider more general case, when the Green 
function has the form 



]B _ Z[p 2 } sU 



8 B A , 2 = 1,2,3,4 (4) 



Here g % [p] is a certain function of the 4 - momentum p that is again diagonal 
in flavor indices. We shall also denote <7s[p] = m[p\. The Green function of 
this form appears, in particular, in the lattice regularization considered in 
Section [6] when Lorentz symmetry is broken due to the regularization. It is 
worth mentioning that (J3J) is not the general form of the Green function in 
the general case when Lorentz symmetry is broken. 

As the matrices g^ip), Zg(p), are all diagonal we further deal with 

them formally as with numbers. 



2.2 When A/3 is the topological invariant? 

In [13] the following expression has been considered: 

^3 = TT^Tr 7 V / Q~ x dG A dg- 1 A dQ (5) 

Here integration is performed along the hypersurface u = of momentum 
space. This is the topological invariant protected by symmetry, i.e. it is 
invariant under perturbations only when the Green's function at to = (or 
the effective Hamiltonian) commutes with the matrix 7 5 7 4 (for topology pro- 
tected by symmetry see Refs. [91 13] and recent paper [H]). For the nonin- 
teracting fermions with positive masses we have A/3 equal to the number of 
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flavors of Dirac fermions. If the interactions are introduced, the Green func- 
tion is changed: Q — > Q + 5Q. If the interactions do not violate the condition 
[0(0, p), 7 5 7 4 ] = 0, then (JSJ) transforms as follows: 



s & = — 2 / tr(^ 4 {[5g]dg- 1 + gd[5g- 1 ]}Agdg- 1 Agdg- 1 

= f tr {^{-G[8G- x ]QdQ- x + gdisg- 1 ]} A GdG~ x A gdCT 1 

-2^2 /_ o tr {f-f{[8G- x ][dS\ + dt^- 1 ]^} A d^T 1 A dS) 
= --^ j[ =o dtr( 7 5 7*{[^- 1 ]^}d&- 1 A«») 

= "2^a / as tr (tVW" 1 ]^- 1 A dg) (6) 



Here <9£ is the boundary of the hypersurface u = in momentum space. 
That's why A3 is the topological invariant if 

J tr ( 7 5 7 4 [8Q- X \Q dg- 1 A dg^j = (7) 

When the Green function is of the form (J3j), we consider the surface <9£ such 
that p 2 = const — > 00 and denote n J = p l /|p|). Then Eq. ([7]) is equivalent 
to 

8ra , , . 

2 ► 0, p -> 00 (8) 

|p| (1 + ^)2 ' ^ 

Thus, we came to the following 

Lemma 

In i/ie 4.D system expression |3]) zs £/ie topological invariant if [Q (0, p) , 7 5 7 4 ] = 
and either momentum space is closed or on its boundary Eq. Fty holds. If 
the Green function is of the form of Eq. |3|) then the latter condition is 
reduced to Eq. fiEj). 



2.3 The asymptotic free models 

It is natural to suppose, that in the realistic theory in deep ultraviolet at 
p 2 — > 00 the system becomes effectively massless and chiral symmetric. We 
imply that the model is assymptotic free and the fermions are coupled to the 
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nonabelian gauge field. For a recent numerical investigation of the fermion 
propagator in such models see, for example, [TS]. Numerical investigations 
confirm that m ^ p < — > while Z{p 2 ) — > at p 2 — > oo. Therefore, (JHJ) 

is fulfiled. Then &A/3 = and A/3 = Sp 1 = K (the trace here is over 
the flavor indices) is equal to the number flavors of Dirac fermions until the 
phase transition is encountered. It is worth mentioning that in the case when 
Z{p 2 ) — > and m{p 2 ) — > const p 2 at p 2 — > 00 we also have 5J\f 3 = (see 



3 Topological invariants for hD massless fermions 

3.1 Topological invariant in 5D at u = 

We consider the system with massless fermions. Let us introduce the follow- 
ing expression: 

W z = ^Tr 7 4 / Q~ x dQ KdQ~ x NdQ (9) 

Here integral is over the 3D hypersurface that encloses the Fermi point in 
momentum space in the u = subspace. Expression is invariant under 
small perturbations of the Green function if {5Q(0, p),7 4 } = 0. That's why, 
the Fermi point in 5D is protected by the given topological invariant if the 
interactions are such that {5Q(0, £>),7 4 } = 0. For example, if the Green 
function is of the form: 



g(u, P ) = ,,„ 4 ^ ; - , a = 1,2, 3, 5 (10) 



z(f) 

^7 4 + PaJ 

then we obtain M% = 2 Sp 1 (the number of 5D Weil fermions). 



3.2 Topological invariant in 5D space (w,p) 

Below we generalize the construction of the topological invariant J\f% consid- 
ered above. The resulting construction uses the Green function defined on 
the surface that encloses the Fermi point in oj — p space. Suppose that the 
5.D Green function has the form 

0( w ,p) = ^i a = 1,2, 3, 4, 5, p 4 = u (11) 

9a\P\T 
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with real functions g a [p\. 

Let us define the function in momentum space 



n = , (12) 

Here the trace is over the Dirac indices only. We can express % as % = 
9al a , 9a = \g\ = y/g a g a , a = 1,2,3,4,5. Now let us consider the 

following integral over the closed 4D surface £ in momentum space such 
that Q does not have poles on E: 



Ma = -4r^ Tr / UdUAdUAdUAdU 
16^4! Jt, 

' J> :e abcde [ g a dg b Adg c Adg d Adg e (13) 



4vr 2 4! 



The given expression for the invariant is reduced to Q for the Green 
function of the form Without interactions Q has the pole at p = to = 
that corresponds to the Fermi point. In this case Ma = 2 Sp 1 for any surface 
that encloses the pole. 



4 Dimensional reduction of the 5D construc- 
tions 

4.1 Dimensional reduction of N$ 

Let us consider the 4D system with the Green function of the form (T4|). We 
have 

2 = ■ ,r i ff} r i 5 = ( 14 ) 
l 9 \P\1 7 + rn[p\'~) 

Next, we define new 7 - matrices T 5 = 7 5 , T l = i^j 5 that satisfy {r a , T b } = 
25 ab , a = 1,2,3,4,5. We also denote gs\p] = m\p]. Then Q has the form 

2=^t' a = 1,2,3,4,5 (15) 

Now we can treat Q as a 5D Green function defined on the hypersurface 
P5 — 95\Pi,P2,P3,Pa\- We may consider, for example, the following analogue 
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of expression Q: 

A/*3 = -^TrT 4 / Q-'dQAdQ- 1 AdQ 
6 24tt 2 A,=o 

= 2^ Tl ' l5 ' li L-o g ~ ldgAdg ~ lAdg = ^ 3 (16) 

Here integral is over the closed 3D hypersurface in the slice u = of 4D 
momentum space. One can see that we come to the invariant ^3 discussed 
in Section [2j This equivalence gives the way to calculate A/3. For example, 
let us consider the case of noninteracting massive fermions. We consider the 
closed surface £ composed of the two pieces Ei : p 5 = m and £ 2 : p$ = —m 
connected at infinity by the additional piece of surface. This surface encloses 
the pole of Q in 5D momentum space. Then we calculate expression (Q 
for this surface N% = 2Spl. Therefore, A/^m) — Nz{—m) = 2 Sp 1 while 
Afzim) = —Nz(—m). As a result Nz{m) = Spsignm. 



4.2 Dimensional reduction of A/4 

Let us define the function in momentum space 



M "Tfe^ r ° (17) 



(Remind that here g a is the matrix diagonal in flavor indices.) 

Now let us consider the following integral over 4D momentum space: 



Afi = ^r-rTr / n' dH' A dW A dW A dW 
167r" ! 4! J 

= e abode J 9 a dg b A dg c A dg d A dg e (18) 



If momentum space Ai is closed (for example, has the form of torus or 
4D - sphere), then expression (fl8l) is the topological invariant that measures 
the degree of mapping g : A4 — > S 4 ® S 4 ... <8> S 4 . For the calculation of A/4 in 
open momentum space see discussion below in Section 15.11 



4.3 Dimensional reduction of A/5 

Now let us consider the other way to derive the topological invariant for the 
4D gapped systems starting from a 5D construction. Let us consider the 
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Euclidean Green's function on the 4D lattice Q as the inverse Hamiltonian 
in 4D momentum space and introduce the 5D Green's function: 

G _1 (p 5 ,P4,p) =P57 5 + ^ 1 (P4,P) = {ips + Q _1 (P4,p))(-«7 5 ) ■ (19) 
Then one can introduce the topological invariant as the 5-form: 

N 5 = ^-^Tr / GdG- 1 A GdG- 1 A GdG' 1 A GdG' 1 A GdG' 1 , (20) 
2ii d 5h J 

where the integration is over the Brillouin zone in 4D momentum space (p^, p) 
and over the whole pa axis. In this form the invariant is applicable also to 
the interacting case. Similar to Section |2] we obtain that expression ( 1201) is 
the topological invariant if the variation of the Green function 5Q satisfies 

/ trim-^GdG- 1 AdG AdG- 1 AdG) =0, pi ->• oo (21) 

Jd[M®R] V ' ' 

Remarkably, here Q may be almost arbitrary. In particular, if M. is compact 
closed space, then (12T|) is satisfied for any Q that does not have poles. If M. 
is noncompact, then (l2~Tj) is satisfied for the Green function of the form (141) 
when the system is massless in ultraviolet. 

In terms of the 4D Green function expression f)20p is 

1 r°° r 1 1 

= VTJr / rf ^Tr / T 5 ^^ dQ' 1 

2-K-Hh J-oo Jm P5T + y 1 P57 s + £ 

P5T + y Pbi + y Pbi + y 

In the case, when the 4D Green's function has the form of Eq. (J5J): 
Q = z j p] . , we denote Q = -j^, b = 1,2, 3, 4, 5, and f a = £. Then 

Ad/r a (/ 2 + pi) + (/r - i P5 )d/ 2 ) a dfr (23) 

The integral over p 5 is equal to the residue of the corresponding pole at 
P5 — if that can easily be calculated. As a result invariant (j23p after the 
integration over p^, is reduced to the degree of mapping of the 4D Brillouin 
zone to the 4D sphere of unit vector: 

A/" 5 = ^e abcde J g a dg b A dg c A dg d A dg e = M' A (24) 
We proved here the following 
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Theorem 

Eq. defines the topological invariant for the gapped AD system with 

momentum space A4 if Eq. I[21\) holds. This requirement is satisfied, in 
particular, for the system with compact closed M. and for the system with 
the Green function of the form |3|) that is massless in ultraviolet. 

For the system with the Green function of the form expression for the 
topological invariant Eg. is reduced to Eq. §2ljty . 

The 5-form topological invariant fl20|) has been discussed in [9J [23j 121] . 
In particular it is responsible for the topological stability of the 3+1 chiral 
fermions emerging in the core of the domain wall separating topologically 
different vacua in 4+1 systems (see Sec. 22.2.4 in [9]) and fermion zero 
modes on vortices and strings |21j. The topological invariant for the general 
2n + l insulating relativistic vacua and the bound chiral fermion zero modes 
emerging there have been considered in [TJ El EE2] • 

4.4 Topological invariant A/4 

In addition to the invariant A/5 let us also consider a different construction 
that coincides with A/5 for the case of free fermions 

A? 4 = — *— Tr 7 5 / dg- 1 AdQA dg- 1 A dQ (25) 
4Stt z Jm 

Here the integration is over the whole 4D space J\4. The expression in this 
integral is the full derivative. Therefore, we have: 

a? 4 = -^Tr 7 5 / dg- 1 a dg a dg- 1 a dg 

4Stc z Jm 

= -^Tr 7 5 f g-Hg a dg- 1 a dg 

4871^ JdM 

Here the integral is over the 3D hypersurface dAi. The last equation is 
identical (up to the factor 1/2) to that of for the invariant A3 for massless 
fermions. Therefore, Eq. (|25|) defines the topological invariant if the Green 
function anticommutes with 7 s on the boundary of momentum space. 
For the noninteracting fermions 

£(^,p)^^-,p 2 ^oo (27) 

pa 
10 



(26) 



Therefore, A/4 = Sp 1 (the number of Dirac fermions). The value of the given 
invariant coincides with A/5 in this case. However, it will be shown in Section 
[6] that in the other nontrivial cases like the lattice regularized models the 
two invariants give different values. In particular, A/4 = while A/5 may be 
different from zero for the systems with closed momentum space A4. 

When the interactions are turned on, but in the deep ultraviolet at p 2 — > 
00 the system becomes effectively massless and chiral symmetric, then Q(p) 
anticommutes with 7 5 at p 2 — > 00. As a result, the value of A/4 remains equal 
to the number of Dirac fermions until a phase transition is encountered. 

5 Topological invariant A/5 
5.1 A way to calculate A/5 

In this subsection we consider the Green function of the form (j3J). Let us 
introduce the following parametrization 

g 5 = cos2a, g a = k a sm2a (28) 

k may be undefined at the points of momentum space A4, where g a = 
0, a = 0,2,3,4. In nondegenerate case this occurs on points p^ i = 1, .... 
We denote the small vicinity of the regions where k is undefined by Q = 
^(po) + ^(pi) + •••• Then the expression for A5 can be rewritten as follows: 

3 r 

Af 5 = -^—re abcd / sm 3 2ada A k a dk b A dk c A dk d 
ii 2 A\ JM-n 

= \-^abcd I (3cos2a - cos 3 2a)A; a dk h A dk c A dk d 

n z 4:l JdM-dn 

= -^abcd I (3^5 - 9l)k a dk b A dk c A dk d , k a = -4^9) 

For the convenience we denote all boundary of momentum space (that 
may be consisted of several pieces) by dM. = —dQ^poo). This allows to 
rewrite the last expression as: 

A/5 = ^TT E tabcd I (3g 5 -g!)k a dk b Adk c Adk d (30) 
71-241 , nTl.x -W) 

Let us consider the case when g 5 = on dM.. In particular, if the bound- 
ary of momentum space is placed at infinity and the system is effectively 
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massless in ultraviolet, then g 5 [p] = at \p\ — > oo. We also take into account 
that (75 = ±1 at Pi,i 7^ oo. That's why we obtain: 

^5 = 4 E ^bcd I sign(g 5 )k a dk b Adk c Adk d (31) 



t=0,l, 



Here the sum over p { does not include p^. 

In the case when there is only one zero of g a (we call this point p ) we 
also obtain: 



= I sign (g 5 )k a dk b A dk c A dk d 

7r 2 4! Jdn( Po ) 

= -^r.tabcd I sign(g 5 )k a dk b Adk c Adk d (32) 

7l z 4\ JdM 



Suppose that the effective mass m\p] does not change sign anywhere. 
Then in this case A5 is equal to the degree of mapping 

k a = €M : S 3 (\p\ ^oo)^S 3 ® ... ® S 3 



^[g l {pW + [g*{p)Y + [g 3 { P )Y+[g\p)Y 
Here we proved the following 



(33) 



Theorem 



In the system with the Green function of the form the topological invariant 
A/5 is given by the sum in Eq. I[3l\) if g^,[p\ = on the boundary of momentum 
space. 

When the Lorentz symmetry is not broken at \p\ — > 00, and g a [p] ~ p a , 
we have A/5 = Spsignm, where m is the mass. That's why for the free 
fermions with positive masses A/5 is equal to the number of the flavors of 
Dirac fermions. 

The situation, when several points pi appear within the Brillouine zone 
will be considered later when we shall calculate A/5 for the lattice fermions 
in Section [6j 

It is worth mentioning that the method suggested above could be easily 
extended to the calculation of A3. 
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5.2 Index theorem 



In this subsection we prove the following 

Theorem 

Suppose that the 4D system depends on parameter (3 and there is a phase 
transition at (3 C with changing 0/A/5. We denote A/5 = n + for (3 > (3 C and 
A/5 = n_ for (3 < (3 C . In addition we require that Q does not have zeros and 
does not contain 75 matrix (i.e. Tr(? _1 7 5 = 0). Momentum space Ai of the 
4D model is supposed to be either compact and closed or open. In the latter 
case we need that Q does not depend on (3 on dM. . Then at (3 = (3 C there are 
nf = (n + — n_)/2 flavors of massless Dirac fermions. 

Let us consider 6D space (pi,P2,P3,P4,P5,Pe) with p 6 = (3. 4D subspace 
P\iP2iPzi Pi has the form of the 4D Brillouin zone M while p§ axis and p§ axis 
are straight lines. We then consider the 5D hypersurfaces S + at p§ = (3 + > (3 C 
and £_ at p& = (3- < (3 C . These surfaces, in turn, contain the Brilloiun zones 
of the initial model at ps = for (3 = (3 + ,(3- respectively. We connect £± 
at p 5 = ±00 and at the boundary of M. (this is possible because Q does not 
depend on (3 at dM). The resulting hypersurface £ is closed in 6D space. 

Let us construct the QD Green function 

G~ 1 (P6,P5,P4,P) = Psl 5 + G~ 1 (/3,P4,p), P6 = (3 (34) 
Then we consider the invariant 

A4| E = — |-^Tr / GdG- 1 A GdG- 1 A GdG' 1 A GdG' 1 A GdG' 1 , (35) 
27r d 5!« jy, 

We may neglect in this integral the region placed at infinite p$. That's why 
A/" 5 |s = Af 5 (f3 + ) - A/K/3-) = n + - n_ (36) 

On the other hand, we may deform S in such a way that it is not already 
placed at infinity. A/s|s = n + — n_ means that the 6D system described by 
Green function ( |34l) has (n+ — n-)/2 massless Dirac fermions. 

Q does not contain 75 matrix. Therefore, there are no poles of G at 
nonzero p§. That's why we may deform S in such a way that on this hy- 
persurface £>5 = ±e with arbitrary small e. Also we make f3± infinitely close 
to (3 C . As a result the 6D model with Green function (|34"|) and (n + — nS)/2 
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flavors of massless Dirac fermions is dimensionally reduced to the AD system 
with (3 = f3 c . There are (n+ — n_)/2 poles in (I3"4"j) at p G = (3 c ,p 5 = 0. We 
conclude, therefore, that there are {n + — n_)/2 poles in {?(/3 c ,P4, p)- 

The example of the AD system with compact momentum space that satis- 
fies the given theorem is given by the model with Wilson fermions considered 
in Section O The example of the system with noncompact momentum space 
is given by the free continuum fermions (in this case mass plays the role of 
parameter f3). 

It is worth mentioning that a similar theorem can be proved that relates 
the jump in A3 with the number of massless fermions that exist in the u = 
hypersurface of the Brillouin zone within the system at the interface between 
two gapped states. (Illustration of this theorem see also in Section E]). 

6 Momentum space topology of lattice regu- 
larized model 

6.1 A/3 for free Wilson fermions 

In this section we consider the case of the single flavor of lattice Dirac 
fermions. In lattice regularization the periodic boundary conditions are used 
in space direction and antiperiodic boundary conditions are used in the imag- 
inary time direction. The momenta to be considered, therefore, also belong 
to a lattice: 

2ixK a 27rfQ + 7r . . 

p a = — — p 4 = , K a ,K 4 eZ a = 1,2, 3 (37) 

Here N a , N t are the lattice sizes in x, y, z, and imaginary time directions, 
correspondingly. However, for the lattice of infinite volume we recover con- 
tinuous values of momenta that belong to the AD torus. 

Let us consider the simplest regularization of fermion action called Wilson 
fermions. In the absence of interactions the Green function has the form: 



(J2 r a sinp a - z(m + 5^(1 - cosp a )) 

a a 

Ea r a sinp a + i(m + £ a (l - cosp Q )) 

Ea Sm 2 p a + (m + £ (l "COS p a )) 2 ' 



1,2,3,4 



(38) 
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m = Ma 


M 1 


(2/a + M) 4 


(4/a + Mf 


(6/a + Mf 


(8/a + M) 1 


m = -2 + Ma 


(-2/a + M) 1 


M 4 


(2/a + Mf 


(4/a + M) 4 


(6/a + M) 1 


m = -4 + Ma 


(-4/a + M) 1 


(-2/a + M) 4 


M e 


(2/a + M) 4 


(4/a + M) 1 


m = —6 + Ma 


(-6/a + M) 1 


(-4/a + M) 4 


(-2/a + Mf 


M 4 


(2/a + M) 1 


m = -8 + Ma 


(-8/a + M) 1 


(-6/a + M) 4 


(-4/a + M) e 


(-2/a + M) 4 


M 1 



Table 1: The spectrum of the system with free Wilson fermions. In the first 
column the values of m are specified. In the other columns masses of the 
doublers are listed. Everywhere it is implied that \Ma\ « 1, where a is the 
lattice spacing. Expression v x means x states with masses equal to v. 



m 


A/3 


A/5 


m > 








-2 < m < 


-2 


-2 


-4 < m < -2 


4 


6 


-6 < m < -4 


-2 


-6 


-8 < m < -6 





2 


m < -8 









Table 2: The values of topological invariants A/3 and A/5 for free Wilson 
fermions. 

In Table Q] we represent the spectrum of the model for different values 
of m. The diagonal values in this table represent physical massive states 
while the off - diagonal elements of the table represent unphysical doublers 
with masses that tend to infinity in the limit a — > 0. Due to the periodical 
boundary conditions A/3 remains the topological invariant on the torus. In 
Table [2] we represent the values of A/3 versus the values of m. Here A/3 = 
\m+ll\ ' wriere the sum is over the fermion doublers in 3D, m + 2i 
is the mass of the i-th doubler while C3 is its degeneracy (for the derivation 
see [8]). It is worth mentioning that on the given lattice A/4 = 0. The value 
of A/4 may be different from zero if Q has poles or zeros somewhere in AD 
momentum space. 

Let us notice that the values of the considered topological invariants on 
the lattice with periodical boundary conditions contradict to the continuum 
result. This is due to the fermion doublers that give essential contributions 
in spite of their unphysically large masses. 
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6.2 A/5 for free Wilson fermions 

As follows from Table [2j the invariant A3 does not resolve between the fully 
gapped states in the region —8 < m < —6 and in the region m < —8, which 
are separated by the gapless state at m = — 8 and thus should be topologically 
different. The reason for that is that the invariant A3 is determined only at 
u = 0, while at m = —8 the only gapless state has nonzero u (that is cos£>4 = 
— 1). Since at m — —8 there are no gapless fermions with u — 0, there is 
no jump in A/3 across this transition. However, the topological invariant, 
which is determined for all u, may differentiate between the two states with 
A3 = 0. The invariant A5 plays such a role, and one has the index theorem 
(see Section I5"72j) . which states that the total number n F of gapless fermions 
emerging at the critical values of mass m is determined by the jump in A/5. 
While the jump AA3 across the transition determines the number of those 
gapless fermions which live in the vicinity of zero frequency, n F (u! = 0). The 
resulting index theorems read: 

n F = ^AAf 5 , n F (u = 0) = ^AAf 3 . (39) 

According to Section this equation is valid for the interacting systems too. 
In some cases, however, when momentum space has a complicated form, one 
must take into account not only poles of the Green's function of massless 
fermions but also zeroes [161 122] . 

For the lattice model with Wilson fermions we use the method described 
in Section f5 . 1 1 in order to calculate the invariant A/5. One obtains that the 
zeros of g a appear at the positions of the fermion doublers. The corresponding 
values of sign g§ coincide with the signs of the masses of the doublers. The 
positions of doublers are p ni = (irui, 7m 2 , vm 3 , 7m 4 ), rii = 0,1. We also 
have dbg a (p ni ) ~ (— l)""^-^— j- , where m ni is the mass of the doubler. That's 
why, using expression f l3T|) we obtain: 

Ar 5 =Al = E(-im^|^, (40) 

where the sum is over the fermion doublers in 4D, m + 2k is the mass of 
the k-th doubler while C\ is its degeneracy. This is similar to the result 
obtained in [8] for the different space-time dimension. The values of the 
invariant A5 as a function of bare mass m are also represented in Table |2j 
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The index theorem Eg. ( 13 9 p gives at the transition points m = 0, —2, —4, 
—6,-8 the number of gapless fermions correspondingly rip = 1,4,6,4, 1 and 
n F (uj = 0) = 1,3,3,1,0. 



6.3 Turning on interactions 

When the interaction of Wilson fermions with the lattice gauge field IA = e zA 
defined on links is turned on, we have (in coordinate space): 



g(x,y) = % - J DUexp(-S G [U\) Bet(V[U,m])V^ y [U,m] (41) 
where Sq is the gauge field action while 

V x>y [U, m] = - - Y^[(l + f)5 x+e ^ y U x+ei >y + ( 1 - 7') 5 x _ ei ,yU x - ei ,y} + (m + 4) 5 xy 

(42) 

Here ej is the unity vector in the % - th direction. We expect that the 
Green function in momentum space has the form [T9] : 

m (43) 



with unknown functions Z\p], g a [p], and m[p\. That's why the Green function 
for the interacting Wilson fermions has the form flj]) as well as the Green 
function for the noninteracting fermions. 

Free Wilson fermions exist in several phases for different values of m. 
These phases are marked by the values of A/5. When the interactions are 
turned on, the value of A/5 remains equal to its initial number until the phase 
transition is encountered. This is related to the fact that zeros of g a [p] cannot 
disappear without changing the phase structure of the model. These zeros, in 
turn, are related to the value of A/5 according to (|22p . It is worth mentioning 
that in order to consider interactions with nonabelian gauge fields we must 
consider several flavors of the fermions as it was done in the previous sections. 
All results obtained in the present section may be generalized to this case in 
a straightforward way. 
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7 Conclusions and discussion 



In the present paper we discuss the topological invariants in momentum 
space of the 4D models with the mass gap. In particular, the expression 
for A/3 introduced in [13] is considered. We find the conditions under which 
A/3 is the topological invariant. We also introduce new functional A/5 that 
is represented as an integral over the whole 4D momentum space. A/5 is 
obtained via the dimensional reduction of the 5D construction. In the form 
of Eq. ( 122]) this functional appears to be the topological invariant for the 
wide class of models. In particular, it is the invariant for the asymptotic free 
gauge theories and for the models with compact closed momentum space. In 
addition, we prove the index theorem that defines the number of massless 
fermions % in the intermediate vacuum, which exists at the transition line 
between the massive vacua with different values of A/5. 

In order to illustrate the properties of the mentioned constructions we 
consider the lattice model with Wilson fermions. We found that the vacuum 
states of lattice models with fully gapped fermions (insulating vacua) in 4D 
space-time are characterized by two topological invariants, A/3 and A/5. Due 
to the index theorem mentioned above, they are responsible, in particular, 
for the number of gapless fermions which appear at the topological transition 
between the massive states with different topological charges. The vacuum 
in the intermediate state represents the topological semimetal. According to 
Eq.( )39|) the jump in A/5 determines the total number of massless fermions 
in the intermediate state, while the jump in A3 gives the number of those 
massless fermions, which have u = 0. 

The simple pattern described above may have a certain connection to the 
unification of fundamental interactions. Namely, masses of all elementary 
particles are extremely small compared to the Planck scale, which may indi- 
cate that we live in the vicinity of the critical line of quantum phase transition 
like the transition between the two gapped states in the model with Wilson 
fermions. There is a special reason why nature may choose the vicinity of 
the transition line: the massless fermions in the intermediate vacuum states 
are able to accommodate more entropy than the gapped states [13] . 

It would be interesting to apply the analysis similar to the presented 
here to the consideration of the other lattice models with more complicated 
structure of momentum space (for example, to the overlap fermions [26J). 

Also it would be interesting to consider lattice models that describe vacua 
with broken symmetry, such as superconducting vacuum in Ref. [22], tech- 
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nicolor models with spontaneously broken chiral symmetry etc. However, 
such an analysis may necessarily involve higher Green functions because the 
two - point functions only do not describe appearance of the condensates 
completely 
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